Self-organization, Localization of Shear Bands and Aging in Loose Granular Materials 
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We introduce a mesoscopic model for the formation and evolution of shear bands in loose gran- 
ular media. Numerical simulations reveal that the system undergoes a non-trivial self-organization 
process which is governed by the motion of the shear band and the consequent restructuring of 
the material along it. High density regions are built up, progressively confining the shear bands in 
localized regions. This results in an inhomogeneous aging of the material with a very slow increase 
in the mean density, displaying an unusual glassy like system-size dependence. 



PACS numbers: 45.70.-n, 45.70.Mg, 05.65. +b 

A large class of materials are handled in the form of 
dispersed solid grains at some stage of their processing. 
Thus the description of the rheological properties of sus- 
pensions, pastes and dry granular media is a key question 
which controls the ability of mixing, storing, transport- 
ing etc, these disperse media (l]-^]. Granular systems 
constitute an intermediate state of matter between flu- 
ids and solids 0,^|: they flow like fluids but they also 
build piles indicating that a non-vanishing static shear 
stress is present which is characteristic of solids. From 
this point of view it is also of major interest to under- 
stand the shearing process in these systems. A number of 
experiments have been carried out on the shear process 
in granular materials HQ. Most of these are triaxial 
tests |?],H to determine macroscopic properties such as 
the shear stress or the volumetric strain, as a function of 
the shear strain. 

The intimate interplay between the geometrical ar- 
rangements and the frictional properties of the grains de- 
termines the precise form of the rheological behavior to 
be used at a continuum level. The underlying question 
is the identification of relevant internal variables. The 
most obvious one is the density of the sample, which can 
be made to vary over a wide range by the method of 
preparation. Compared to other parameters describing 
the texture (e.g. fabric tensors accounting for the distri- 
bution of contact orientations) the density has the most 
drastic impact on the stress needed to shear the mate- 
rial as well as on the mode of shearing; from an appar- 
ent homogeneous strain for loose packings to a localized 
steady shear band for dense assemblies || . The coupling 
of the density to the shear properties can be understood 
through the concept of dilatancy Q . 

A related question is whether statistical fluctuations 
have an impact on macroscopic properties. Lately, there 
has been an upsurge of interest in trying to characterize 
the large stress fluctuations pC|-[T^] in silos, Couette flow 
or slider block geometries, or to understand the statistics 



of interparticle contact forces Recently, spectacu- 

lar experiments in two-dimensional Couette shear cells 
were carried out ]l3| where the movement and stress of 
individual particles were monitored in order to describe 
the inner structure and the force network in the sheared 
granular material. It was demonstrated that stationary 
motion is accompanied by large stress fluctuations due 
to the formation and breakdown of arches. Large fluctu- 
ations were also found in three dimensional steady state 
shear cells [ p"5[ . 

This issue has also been raised by the results of re- 
cent numerical simulations of rigid grain assemblies [|l6| , 
where even at low densities, the shearing which appears 
as homogeneous over long times, in fact consists of a 
succession of sudden changes of quasi-instantaneous and 
localized strain fields. This observation suggests that the 
transition from the particle based description to the con- 
tinuum one requires the detailed understanding of the 
statistical features associated with these sudden changes. 

In this Letter we present a simple model for the shear- 
ing of a granular medium in loose samples. We describe 
the strain field at every instant as a shear band, chosen 
through a global optimization procedure, which is equiv- 
alent, as we shall see later, to searching for the ground 
state of a directed polymer in a random potential Jl7],|l8| . 
However, this potential is not a priori frozen in but has 
a self-organized development due to our procedure of 
choosing and changing the shear band. Though very sim- 
ple and with only the minimum of ingredients, the model 
shows that the density of the medium increases anoma- 
lously slowly. Further we are also able to predict on the 
basis of this model, that large scale inhomogeneities build 
up in a system subject to a steady shear. This could be 
an interesting feature to compare with experiments. 

Let us consider a shear process, assumed to be invari- 
ant along the shear direction [z in Fig. ^). This geom- 
etry is appropriate for instance, in an annular shear cell 



of large radius |15|. We consider moreover a continuum 
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FIG. 1. Schematic picture of the shear process. The shear band 
is parallel to the shear direction z due to periodic boundary condi- 
tions in this direction. 

description, valid on scales much larger than that of 
individual grains. We now introduce a fundamental as- 
sumption of our model: We assume that the instanta- 
neous strain field is always localized on a single shear 
band j|,[l9| . Experimentally, it is known that shear bands 
have a typical width of about ten grain diameters. Thus 
at a continuum level, the velocity field is indeed discon- 
tinuous across the shear band. From the geometry of 
our set-up, the shear band must be a continuous surface 
due to topological constraints (Fig. [l]). Further, we as- 
sume that because of the translational invariance along 
the z axis, the system can be reduced to a two dimen- 
sional one in the x-y (cross-section) plane, through an 
averaging over the z direction. 

The basic hypothesis of the localization of the shear on 
the shear band at all times, is not as restrictive as it may 
appear. We only refer here to instantaneous shear rates, 
and provided the shear band changes rapidly enough, 
coarse-graining the strain field in time may produce a 
uniform shear rate. Experimentally, though it is very 
difficult to have direct access to the instantaneous shear 
rate, large fluctuations found in the shear stress may in- 
dicate that the shear is never quite uniform, even at early 
times. As mentioned earlier, this seems indicated also by 
numerics p6| . 

Initially we consider a loose-packed sample. At a suit- 
ably coarse-grained scale the medium can be described 
as a continuum, where the density is a random function 
displaying fluctuations around a mean value. Under a 
constant normal load, a threshold shear force (or torque 
for an annular shear cell) has to be applied to impose 
a non-zero strain. Locally, after integration along the z 
axis, the density controls the threshold shear force. Al- 
though this is inessential, for simplicity we assume that 
the ratio of shear to normal stress, i.e. the friction co- 
efficient, increases linearly with density. As mentioned 
earlier, the texture of the medium also contributes to the 
friction coefficient. However, since we consider only shear 
in a fixed orientation, a single scalar parameter combin- 
ing density and texture should suffice. This parameter is 
called "density" for short and is denoted by g(x, y). Thus 
at any time the state of the medium is characterized by 
this field. 



We determine the shear band (path in the (x, y) plane) 
by the following three conditions: a) it is continuous, b) 
it spans the sample in the x direction without overhangs 
and c ) the sum of the density along it is minimal among 
all possible paths satisfying a) and b). One can recognize 
that this is the well known problem of finding the ground 
state of a directed polymer in a random potential [ljj ] . 

Relative motion of the particles takes place within the 
shear band while the rest of the sample remains still. 
Small movements can totally rearrange the local struc- 
ture Jl5|,|2(| and thus may induce large changes in the 
local density. We simplify this complex behavior by re- 
newing the density g only along the shear band, by inde- 
pendent random values taken from a fixed distribution. 
After this, a new shear band is again located as described 
above. Thus the shear process consists of a succession of 
localized slips occurring at very small time scales. We 
note that in characterizing this process, in the spirit of a 
continuum modeling, we ignore potential stress inhomo- 
geneities in the medium. It is a simplifying assumption 
of the model to relate the shear band localization only to 
the density, and not to the full solution of the local stress 
distribution. 

In order to be able to simulate the above model we dis- 
cretized it on a square lattice either with principle axis 
parallel to x and y and considering first and second near- 
est neighbours, or tilted by 45° considering only nearest 
neighbours. Periodic boundary conditions are imposed in 
the y direction. Simulations with site and bond versions 
were also carried out leading essentially to the same re- 
sults. We consider here square samples with system size 
N x N with N varying from 32 to 512. Initially a den- 
sity Qi (a random number uniformly distributed between 
and 1) is assigned to every bond i. We define the in- 
stantaneous shear band as the spanning directed path 
along which ^2 Qi is minimal ( applying the usual trans- 
fer matrix method JnJ). Once the shear band is found 
the bonds belonging to it are assigned new values taken 
from the same uniform distribution as used initially. We 
repeat this process and monitor different properties of 
the system 

We define the average density (g) as the mean value of 
the density of the sites not belonging to the shear band. 
This definition, as well as our procedure of choosing the 
least and changing it, guarantees that the average density 
is a monotonically increasing function of time. 

The monotonic behavior and the bounded nature of 
the average density (g < 1) ensure that it has an asymp- 
totic value. In finite samples this is equal to 1. In Fig. 
U we have plotted the deviation of the average density 
from this asymptotic value. At early times (t/N < 2) the 
rescaled curves go together independently of the system 
size; later non-trivial system size effects can be observed. 
The relaxation to the asymptotic value gets slower as the 
system size increases. 
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FIG. 2. The difference of the average density (g) from its 
asymptotic value is plotted as a function of time t rescaled by the 
system size N. The system sizes are 32, 64, 128, 256 and 512 from 
to bottom to top respectively. 

Since the system evolves entirely through the process 
of choosing and changing the shear band, we have mon- 
itored the following two important quantities related to 
the shear band: The Hamming distance d (which is the 
number of different sites between successive shear bands) 
(Fig. ||a) and the average density of the sites along the 
shear band (qsb) before change (Fig. ||b). It is appar- 
ent from the figure that there is a characteristic time of 
t c i ~ N, below which the distance is essentially constant 
and equal to the system size and the density of the shear 
band is roughly constant. This can be understood quali- 
tatively from the following considerations. Since the very 
first shear band is equivalent to the ground state confor- 
mation of a directed polymer in a random potential, we 
know from this analogy that the mean density along this 
shear band is much less than 0.5 |17|. Once the path 
is refreshed, its mean density increases to 0.5. The next 
shear band tends to be repelled by the previous one since 
there still exist many spanning paths with a lower den- 
sity. Thus at early times, two successive shear bands 
differ completely (Fig. |^a) and the density of the shear 
band remains more or less the same (Fig. |3|b). This ini- 
tial phase should last until on average all sites have been 
refreshed a few times, a number of time steps of the order 
of N. 

The absence of overlap between successive shear bands 
in this early time regime reflects the fact that no well de- 
fined shear band can be observed in loose granular me- 
dia. Experimentally this is connected to the difficulty in 
quantifying fluctuations, when the mean shear strain is 
of small magnitude. So what is observed is seemingly a 
homogeneous shear. 

There is a transition regime up to t c2 — 20iV where 
we still have a good quality data collapse. In this regime 
both curves d and r = 0.5 — qsb start to fall off. The 
decreasing distance indicates an increasing persistence of 
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FIG. 3. a) Log-log plot of the time dependence of the distance 
d for system sizes 32 to 512 scaled together. Both the distance d 
and the time t scale with the system size N. b) Log-log plot of the 
deviation of the density of the shear band before the update from 
its asymptotic value as a function of time for system sizes 32 to 512 
scaled together. 

the shear band. As the average density of the system 
increases (Fig. ||) the density of the minimal path also 
grows and thus the repulsive interaction between two con- 
secutive shear bands progressively fades away. Finally, 
by the end of the transition regime, the interaction be- 
comes attractive and a much slower relaxation process 
takes place. 

The above measurements point to a localization of the 
shear band, induced by the imposed dynamics. In order 
to understand better how this comes about we present 
density snap shots of the system at four different in- 
stances (Fig. H) varying from t/N ~ 4 to 4000. We 
observe that initially (Fig. ||a) the density appears homo- 
geneously distributed. Then progressively high density 
regions become apparent. The mechanism for the forma- 
tion of these regions is the following: As the average den- 
sity increases, the interaction between successive shear 
bands becomes attractive and the path gets restricted in 
space. Small fluctuations of the shear band then lead to a 
density increase in this region. The presence of these sur- 
rounding areas of high densities increases the attraction 
of successive shear bands, thus leading to a positive feed- 
back process resulting in regions of finite width and very 
high density where the shear band is trapped in the mid- 
dle, in a "canyon-like" structure (black lines surrounded 
by white in Figs. ||c and d). 

The escape from the above described trap is only pos- 
sible via a jump to another local minimum. The proba- 
bility of such a jump decreases faster than exponentially 
with increasing density. Thus as time grows, the aver- 
age jump size decreases even though large regions with 
relatively small densities remain. The progressive self- 
quenching of the shear band in the system is responsible 
for the anomalous slow increase in the average density. 
This inhomogeneous aging and extremely slow dynamics 
is reminiscent of a glassy behavior. 

In order to get some more insight into the slow dy- 
namics of the system we have studied the same model on 
a hierarchical lattice. The simple geometry allows for a 
detailed analytic treatment of the model. This study will 
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FIG. 4. Snapshots of densities at different time on a system 
of size 256 by 256: a) t = 10 3 ~ 47V, b) t = 10 4 ~ 40JV, c) 
t = 10 5 ~ 4007V, d) t = 10 6 ~ 4000/V. The greyscale is indicated 
at the bottom of the figure. 

be reported elsewhere p^] . Here we only summarize 
the main features of this analysis. The slow density in- 
crease and strong system size dependence seen on the 
square lattice are also seen in the hierarchical one. Here 
we can show that 1 — (g) decreases as a sum of power-laws 
with a vanishing exponent depending on the lattice size, 
i.e. the number of generations of the hierarchical lattice. 
Further, the early time regime is a single function of t/N 
as for the square lattice, while the late time regime scales 
instead as t/N a , where a = 1/ log(2). 

In spite of its simplicity, the model we have introduced 
displays some interesting consequences of collective orga- 
nization of density fluctuations in a granular assembly. 
Although only time-independent rules are introduced, 
the simulations reveal a slow densification which occurs 
together with a non-trivial patterning of the density in 
the sample. Simultaneously, the shear strain is localized 
on shear bands which acquire progressively a longer and 
longer persistence. The occurrence of high density re- 
gions confining the shear band is a feature which should 
be observable using X-ray tomography as recently per- 
formed in triaxial tests by Desrues et al || . 
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